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Abstract. We propose a new method to compute the configurational entropy of glassy systems as a function 
of the free energy of valleys at a given temperature, in the framework of the Stillinger and Weber approach. 
In this method, which we call free-energy inherent structures (FEIS) approach, valleys are represented by 
inherent structures that are statistically grouped according to their free-energy rather than the energy as is 
commonly done in the standard procedure. The FEIS method provides a further step toward a description 
of the relaxational behavior of glassy systems in terms of a free energy measure. It can be used to determine 
the character of the glass transition as well as the mode coupling and the Kauzmann temperatures. We 
illustrate the usefulness of the method by applying it to simple models of glasses and spin glasses. 

PACS. 64.70.Pf Glass transitions - 75.10. Nr Spin-glass and other random models - 61.20. Gy Theory and 
models of liquid structure 

1 Introduction stable states [Q. In particular, the free energy landscape is 

characterized by a large number of valleys which in spin- 



Spin glasses and structural glasses are characterized by 
the presence of a complex structure of stable and meta- 
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glass theory are also commonly referred to as basins, pure 
states, phases or ergodic components. These valleys con- 
tain a large number of configurations (typically growing 
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exponentially with the size of the system) with macro- subtle issue. In mean-field theory, where ergodicity break- 
scopic properties (for example the energy or the magneti- ing holds, the decomposition of the phase space into val- 
zation) that in general depend on the specific valley. The leys has a clear physical meaning. In that case one can de- 
complexity of the configurational space topology reflects fine as configurations belonging to the same valley those 
itself into a complex dynamical behavior commonly re- such that starting from any of them, any other configura- 
ferred to as glassy behavior. These ideas trace back to tion in the same valley can be reached (or get close as we 
more that 30 years ago to a seminal paper of Goldstein |2j, like for continuous systems) in a finite time. In the case of 
who suggested that the dynamics of a supercooled liquid short ranged systems, however, the extension of this defi- 
can be understood in terms of a diffusive process between nition remains quite unclear. In fact, in the case of finite 
different valleys. At low temperature the dynamics slows dimensional systems, even in presence of Replica Symme- 
down since the system gets trapped for long time in a val- try Breaking (RSB), activated processes must be included 
ley and the slow long time relaxation is only governed by in the physical picture. 

the inter-valley motion. An alternative point of view rests on a dynamical def- 

inition of the concept of valley Metastable states or 

Along the same lines, but in a different context, Still- 
valleys V can be defined by introducing the observation 

inger and Weber (hereafter referred as SW) proposed in 

time r and looking at the cumulative probability P v (t) 

the early eighties to identify each valley with its mini- 

that the system starting in any configuration contained in 

mum, called Inherent Structure (IS) and to build up an 

the valley V at time t — escapes from V at or before 

IS-based analysis [3]. Each valley contains all those con- 
time t. The valley is well defined |5| if there exists a rea- 

figurations that map into the IS by a steepest descent 

sonably small significance level p such that the following 

dynamics. On the one hand, such a coarse-grained de- 
inequality is satisfied: 

scription of the phase space may be useful to grasp the 

universal and generic features of glassy dynamics. On the p v ( T ) < ^ (1) 

other hand, replacing each valley with a single configura- 

The escape probability, and hence the partition, in gen- 

tion implies of course a strong reduction of information. 

eral depends on the external parameters such as fields or 

A relevant question then is to understand under which 

temperature and on the system size N. If for a valley it 

conditions the resulting IS-based dynamics is equally well 

happens that 

representative of the off equilibrium behavior of the sys- 

lim lim P v {t) = (2) 
tern. From purely theoretical grounds, the decomposition r^oow^oo 

of the phase space into non-overlapping valleys is always then the valley is absolutely confined [5] in the thermo- 
possible. However, its operative definition may be a rather dynamic limit and corresponds to an ergodic component. 
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This is the typical case for mean-field models since free- the validity of the SW scheme which has the advantage of 

energy barriers between valleys growth with the system being quite easy to be implemented. 

size. The paper is organized as follows. In Section|3we dis- 
For non-mean-field systems valleys are not, in general, cuss the P hase s P ace partitioning and identify the relevant 
absolutely confined and their identification is not trivial, physical quantities. The method to compute the FEIS is 
A possible strategy to deal with those systems is that of described in Section H and is applied to two representa- 
looking for possible simple limiting cases and using them tiv e models, the Sherrington-Kirkpatrick (SK) model and 
to partition the phase space. The SW scheme belongs to the Random Orthogonal (ROM), in Sections Hand El The 
this class. Indeed since at T = any barrier cannot be choice of these model is motivated by their different phase 
surmounted, it is clear that all valleys are absolutely con- space organization. Finally Section © contains conclusions 
fined at T = 0, making their identification simpler. The and some perspectives. 
SW scheme identifies valley with a T — dynamics, there- 
fore we can regard the SW scheme as a T = partitioning. _ _ . . , . . ,, 

2 Partitioning the phase space into valleys: 

In the original SW scheme different valleys are clas- 

the free-energy landscape 

sified according to the energy of the minimum (IS). This 

picture can be somewhat misleading since it plays down In genera i 7 given a decomposition of the phase space into 
the role of entropy. In this paper we shall discuss the non-overlapping valleys, we can write the partition func- 
more physical classification in terms of Free Energy Inker- ^j on g as 
ent Structure (FEIS). We also introduce a direct method, 

Z{T) = V exp(-/3W(C)) = expf-^yCT)) , (3) 
which to our knowledge is new, of computing the free en- ^ ^ 

ergy of valleys. We emphasize that the FEIS classification where Q denQtes the gpin configurationS; H is the Ha mil- 
is more general than the usual IS-energy one. The two tonian &nd p _ x ig thfi inverge temperature . jp v(T ) is 
classifications lead to the same information only if differ- the free energy rf ^ yaUey Vj 
ent valleys have similar entropies, which in general may 

not be the case. S ^ l0 ^E «p(-/W(C))). (4) 

r cev 

Strictly speaking, the partitioning depends on temper- Expectation values of observables can be defined in a given 
ature, this is particularly relevant for the dynamical prop- vauey y by restricting the Boltzmann measure to configu- 
erties of the systems. We do not address here this problem ra tions belonging to that valley (restricted ensemble 0): 
since we are interested in understanding the role of entropy 

(A) v = V^(C)exp(-/3H(C)+/3^ v (7 1 ))- (5) 

once a partition is given. For this reason we merely assume 
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Note that of the mean-field theory can be extended to realistic sys- 

.-^ tems 0. A key concept describing this mean field behav- 

(A) = 52{A) v exp(-l3r v (T)+l3r(T)), (6) 

v ior is the configurational entropy (also called complexity) 

where T is the total free energy of the system. After these „ , P . „ . ,. . „ ,, . „ 

o c {j, 1 ) dehned as the number of valleys with a given free 

definitions it is natural to consider the entropy of the val- ., , _ „ , 

energy per site / . 1 he configurational entropy counts the 

ley V, defined as number of valleys accessible to the system in the same 

S\>(T) = f3{{E)\>(T) — T\>{T)). (7) wa V as the usual entropy in Boltzmann theory counts the 

number of configurations with a given energy. Therefore, 

In principle the partitioning of the phase space into val- 

the information that it contains has physical meaning if 

leys allows to extend all thermodynamic relations to more 

an equiprobability hypothesis (equivalent to that of equi- 

general valley-dependent ones. As stated above, the rel- 

librium ensemble theory) holds for spin glasses and glasses 

evant question is whether a given partitioning of phase 

in the off-equilibrium regime [7] . This hypothesis assumes 

space has a physical meaning. 

that valleys with a given free energy / have the same 

Recent results about the dynamics of spin glasses have 

weight and can be treated on the same footing as far as 

shown that a modified form of the fluctuation-dissipation 

the dynamical properties of the system are concerned. The 

theorem (hereafter referred to as FDT) is satisfied in the 

implications of such an assumption for the description of 

aging regime [B]- For structural glasses or one-step RSB 

complex landscapes in glassy systems are very strong. It 

spin glasses these violations are asymptotically described 

provides a flat measure for the dynamics of glassy sys- 

by a single time-scale which can be related to an effec- 

tems similar to the Edwards measure for granular media 

tive temperature. Violations of FDT are tightly related 

nnnn. 

to the spectrum of off-equilibrium fluctuations and the 

dependence on a single effective temperature has been ad- In this framework a description of the off-equilibrium 
vocated as the consequence of the cooperative character dynamics in terms of the complexity has been successfully 
of the dynamics [7j . Although FDT violations have been obtained for structural glass models ^5] within the IS for- 
computed in detail only in mean-field systems there is ev- malism 3 where each valley is identified by its minimum 
idence that they yield a physically valid description of the IS. The physical assumption behind the validity and use- 
free energy fluctuations in the aging state also for short- fulness of this coarse grained description (where configu- 
ranged systems The fact that FDT violations can be rations are grouped into different IS) is that the time-scale 
well rationalized in a mean-field framework 9 (the ef- of inter-valley motion is much larger than the time-scale 
fective temperature is related to the function x(q) in the of intra- valley motion. If typical valleys are locally similar, 
Parisi solution of spin glasses) suggests that some features i.e., they can have different energies but similar entropies, 
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then one may simply count them from the probability dis- crete systems. This method will provide us with the ex- 

tribution of the IS energies and obtain the complexity, plicit dependence of the configurational entropy on both 

apart from an unknown constant. This last assumption temperature and free energy. Moreover it directly gives 

can be relaxed by estimating the entropies of the valleys s c (including also undetermined global constants) with- 

by treating fluctuations with the harmonic approximation out using any approximation for the intra-valley entropy 

11311411511611711%] . or by more realistic approximations (such as the harmonic approximation discussed in the last 

which take anharmonicity into account (unfortunately this but one paragraph) , improving the results obtained with 

procedure is not applicable to discrete systems). It should other methods |18ll4llf)ll6ll7ll8l22| . To our knowledge 

be noted that that such a method for evaluating the com- this is the first time that a method to compute the statis- 

plexity in structural glasses stems out naturally from a re- tical distribution of the free energies of the IS valleys has 

cently proposed thermodynamical approach |19I2()I18| to been proposed and tested, 
the liquid glass transition, within the approximation that 

free energy valleys are well described by the corresponding . . fll , . . 

3 Description of the free-energy inherent 

IS valleys. 

structures (FEIS) method 

We note that the utility of the SW decomposition is 

less evident in the case of Full Replica Symmetry Break- We consider the IS decomposition introduced by Stillinger 

ing (FRSB) spin glasses where the structure of valleys and Weber ( Sw ) El, where each valley V is labeled by the 

is more sensitive to temperature and no clear time scale corresponding minimum (IS), and denote its free energy 

separation exist, making the application of a T — de- P er s ^ e by 

composition to finite T questionable. Nevertheless a recent /is = Ja-^V- (8) 

numerical study 1211 shows that the coarse graining from TT r , n i , -n j , • • 

J *— ' b ° Hereafter small letters will denote intensive quantities. 

configurations to IS preserves the physical properties of r\ ^ ■ ^ t p m\ 

° * r j r i (j ur g 0a j lg compute the complexity o c (j,l ) witn- 

the system. We shall see that a statistical description of , , . . ,, . ., r n 

J out making any assumption on the similarity of valleys. 

the free-energy landscape made in terms of the free energy T , , , , c T JT , cj p™ i. 

OJ * OJ fn the context of Lennard- Jones glasses, Speedy |23] has 

of valleys (rather than their energy) is still useful also for , , j. c , 

J y OJ ' proposed to compute the free energy of a valley by adding 

FRSB spin glasses at finite temperature, and that the cor- TT .,, . ,. , c 

r ° to the Hamiltoman a coupling term between the conhgu- 

responding complexity S c ( f, T) gives a proper description ,. , i , , TO , ,. c .. 

^ ° ^ J XJ ' ' ° f f t ration that labels the IS and the run-time configuration, 

of the glassy phase. , . ^ijcj i a--i 

° J —tq, 1 being an appropriately dehned overlap. A similar 

In the next Section we shall propose a method to nu- approach in which the complexity is evaluated by cou- 

merically compute the valley free energy suitable for dis- pling the run-time configuration to a generic configuration 



6 Barbara Coluzzi et al.: A New Method to Compute the Configurational Entropy in Glassy Systems 

of the valley was considered in the first reference in where g(f,T) is the density of IS with free energy density 

By taking the thermodynamic limit first and extrapolat- / which defines the extensive complexity 
ing later to e — > it is possible to restrict the measure 

S c (f,T)=\og{g(f,T)), (II) 

to those configurations C belonging to the IS valley. This 

procedure is laborious, since it requires extensive Simula- and pis(T) is given by ©. The denominator in JTOj! is the 

tions for different values of the coupling e and for each equilibrium partition function 
different IS. 



Here we propose a different strategy based on the prob- 
abilistic definition [21] of the valley free energy. The dy- 



Z(T) = exp(- Weq(T)). (12) 
Equations I|1U|) - (|12[) l ea d to the relation 



namical evolution of the system defines a probability mea- s c (f, T) = 

io g [Hf,r)] +/?(/ _ /eq(T)) _ (13) 



sure pv over the valleys. In the case of an ergodic dynam- 
ics, and assuming that the observation time is Ttot 3> T eq , 
the statistical weight of the single valley is obtainable di- 
rectly from @ and the Boltzmann a priori equiprobability 
hypothesis |24j : 



N 

The key formula in this expression is the probability V{f , T) 
which can be estimated by computing fig from J^J with 
the following procedure. After having thermalizcd the sys- 
tem at temperature T for ttherm Monte Carlo (MC) steps, 
a steepest descent procedure is used every i run MC steps 



t v , ar , , ... for AT run to identify the IS [2U • The interval t run is chosen 

Pv{T) = = exp {-(3Tv{T) + pT c<i {T)), (9) 

Ttot 

to be of the same order as ttherm to ensure that configura- 
where r v denote the time spent by system in the valley tions obt amed at the end of each period are uncorrected. 
V during the total observation time r tot , and T cq (T) the If thc numbcr of different IS is not too large and N mn is 
equilibrium free energy. Alternatively eq. © can be de- taken lar § e enou S h to ensure that each IS has been visited 
rived from the information theory as the ensemble which a substantial number of times, then p is can be estimated 

minimizes the information to specify one valley when we as 
only know that the system is in contact with a heat bath 

(T) = j^- (14) 

at T {unbiased component ensemble) 0. In either case where Njg is the number of times the given IS has been 
the probability to find at temperature T an IS with free found ^ Njs = Fmm p the IS free energy now 

energy density equal to / is given by reads 

/ls( r) = 4,„ g (|£L) +Wn (15) 

■PiS, T) = (6(f - /b)> = X> (/ - As) Pis (T) 

IS The equilibrium free energy / cq (T) can be computed by 

= g (/,r)exp(-W) 

^ exp (—0Nfis) ' performing a different MC run and integrating the energy 
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of the system from infinite temperature limit down to the This approach relies on the hypothesis that all valleys 

working temperature T: with the same IS energy have the same relevance to the 

, , , , . statistical properties of the system, so that the frequency 

f3f cq (T) = J dp'e cq ({3')-s cq (P = 0). (16) 

of visiting an IS with a given energy only depends on the 

Finally from the value of /is it is easy to construct the 

total number of IS with that energy. If the approximation 

histogram V(f,T) and using eq. compute s c (f,T). _, 1M 

(I2U(I is valid, s c (e) can be evaluated, apart from the un- 

This approach has an important difference with the 

known constant fo(T), directly from (|18|l . As the unknown 

usual SW decomposition. In the standard SW method IS 

constant in s c (e) only depends on T (and not on the en- 

with the same energy eis are assumed to occur with the 

ergy) one can check the correctness of the results from the 

same probability. However, in the FEIS method IS with 

superposition of the appropriately shifted curves obtained 

the same free energy fis are grouped together. Compared 

at different temperatures. Usually it is reasonable to ex- 

to (tTU|) . in the standard SW procedure we have: 

pect that this happens in at least two different situations. 
— eis) Pis (T) First, when the temperature is low enough that only con- 

is 

<? e (e) exp (— f3Nf(e, T)) figurations near the bottom of the IS valleys contribute; 

= £is exp(-/3JV/i S (e,T)) ^ 

second, when the IS-valleys are narrow as in REM-likc 

and therefore __ 

models |26|. As we are going to discuss in detail, this is 

\og[V(e,T)] 

s c (e) ~ j^- h /3(/(e, T) — / oq (T)) (18) not the case in FRSB spin glasses, where the size of valleys 

with 



s c (e) = Yims c (f,T). (19) 



must also be taken into account. In such case the config- 
urational entropy is meaningful only when expressed in 
terms of the free energies of the valleys. 



For models with continuous degrees of freedom s c (e) can 
be obtained directly from O by evaluating f(e, T) with, We conclude this Section by noting that the partition 
e.g., the harmonic approximation. function O can be written as 

For discrete models this decomposition is useful when- 
ever the free energy of valleys with IS energy e has a trivial Z{T) = ^ exp(-/3Nf IS ) = j df g(f, T) exp(-0Nf) 

is 

dependence on e: p 

= J df exp(-0JVsP(/,T)) s (22) 

/(e,T)~e + / (T), (20) 

with / (T) only function of temperature. Note that in this w ith the generalized free energy 
case one has 

s c (f,T)=s c (e + f (T)). (21) <P(f,T) = f-Ts c (f,T). (23) 



8 Barbara Coluzzi et al.: A New Method to Compute the Configurational Entropy in Glassy Systems 

In the large N limit, the equilibrium free energy is given 4 Models and Simulation Details 

by the minimum of <P: 



/ eq (T) = ndn#(/,T) = f*(T) — T s c (f*(T), T), (24) 



We have considered two spin glass models: the Sherrington- 
Kirkpatrick model |27| (SK) and the random orthogonal 
model (ROM) [2E| They are both described by the Hamil- 
tonian 



where we used <|12[) and (|22|) . Under the assumption that 
the typical IS energy e*(T) is similar to the equilibrium 
energy e cq (T) then the complexity s c (f*(T),T) is simply where the spins arc Ising spins (o l = ±1) and the J l} 



H = — Jij<7i<7j, (25) 

l<i<j<N 



given by the difference between the entropy of the system are symmetric (J„ = J Jl ) quenched random couplings, 

at equilibrium s cq (T) and the typical IS entropy s*(T) i n the SK model the J l3 are uncorrelated variables with 

(see Ref. |19|2U|). This allows the numerical computation zer0 mean an( j i/jy variance, while in the ROM are 

of s c (f*(T), T) in systems with continuous degree of free- the elements of a random orthogonal matrix, with J u = 

dom such as structural glass models where s*(T) can be an d 

evaluated as the average entropy of a disordered harmonic J k J. k = AS- (26) 

k 

solid (or also within more realistic approximations). 

We have considered both the SK model with binary cou- 

The formalism developed in 19 20 is appropriate for ,. T , 1 ; nrr , .,, „ . ,. n , T N 

^ ' — 1 — 1 plmgs Jij = il/viV and with Gaussian couplings -r{Jij) = 

glassy systems with a clearcut time-scale separation be- ,„ i* T \-Mi i .r T 2 in\ • ■ • •-, 

° J (ztt/N) ^ exp(—NJ(j/2), obtaining very similar results 

tween the inter-valley and the intra-valley motion, so that 



The reason for this choice of the SK and ROM is 



the valley free energy is a well defined quantity. Our ap- n , - , , ., n , 

based on the fact that they describe the two most rele- 

proach does not suffer from this limitation and we shall r „ . , , , . „ T . . , 

vant scenarios of mean-field glassy systems: tfie SK model 

gives the mean-field theory for disordered and frustrated 

magnets with FRSB, while the ROM describes structural 

glasses in the mode coupling approximation |T|. Recent 

work [221 suggests that the study of finite-sized mean-field 

systems is a useful route to investigate activated processes 

in real systems. A statistical analysis of the IS in the SK 

characterizing the glass transition (e.g. the MCT transi- , , . , . . . . T , . . 

model has been presented m 22 21 . In particular, it was 

tion temperature or the Kauzmann temperature). . . _— . ,, . . ,. , ., . „ , 

shown m |21| that the probability distribution of the over- 

In what follows, we will show the powerfulness of our lap -P(g), i.e. the order parameter that describes the FRSB 

method using simple spin models of glassy systems. transition, can be computed between IS instead of equi- 



show that in (|23(l is the relevant quantity to look at 
also in FRSB model, where no clear time-scale separa- 
tion exists at finite temperature. Moreover it is just the 
behavior of at different temperatures which gives infor- 
mation on the nature of the glass transition allowing to 
predict the values of the different relevant temperatures 
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librium configurations (weighting the IS with their prob- following, to restrict the analysis to a smaller range of free 

ability of occurrence in the simulation pjg(T), which is energy values. 

the weight of the corresponding basin at the considered We have also checked that, as expected these re- 
temperature). In this sense, the coarse-graining from con- suits do not depend on the details of the algorithm used 
figurations to IS seems to preserve the physical properties to determine the IS: we have verified that the greedy algo- 
of the system. For the ROM it has been found in [22| rithm (which follows the steepest path in the phase space) 
that valleys are statistically identical as they have a very and aT = MC dynamics produce the same results. For 
small intra-valley entropy, and the usual SW decomposi- the finite T simulations we have used a MC update with 
tion in terms of the energy of the different IS provides a a simple Glauber dynamical rule. To reduce sample-to- 
fairly good statistical description of the potential energy sample fluctuations for the SK model data have been av- 
landscape and the relaxational dynamics. In the SK model eraged over 10 samples. For the ROM we find very small 
(and, by extension, all models with continuous RSB), con- sample-to-sample fluctuations. Already for systems of size 
trarily to the ROM, valleys can be very different, so the N = 32 different samples almost give the same results. 
IS free energy has important contributions coming from This is most probably due to the high degree of corre- 
intra- valley fluctuations. lation among the couplings imposed by the orthogonal- 
We have considered small volumes, typically N = 32 ity requirement. For this reason, and because simulations 
and 64, to avoid having too many valleys. In fact the num- turn out to be very long, the data for ROM reported in 
ber of stationary points of the energy surface Ms is known the figures have been averaged over only two samples. 



to grow exponentially with the system size: 

Ms ~ exp(aiV) (27) 



5 Results and Discussion 



with a ~ 0.2 for the SK model M and 0.3 for the ROM In Figure[TJwe show the complexity s c (e) for the SK model 

■Km „ r , .~ , . c ., .... as a function of the IS energy as defined in the usual SW 

|30J. Because of that, if we want to satisfy the condition OJ 

MsMun ~ 0(1) in d, TV cannot be too large, otherwise decomposition computed from eq. ©. We have consid- 

, r /AT d TC • rn ■ ,i i j t 4.1, ercd two different averages: 

A/is/jVrun = for those IS inefficiently sampled. In the case ° 

of the SK model for N = 32 we have Ms 10 3 while for _ thc annealed ave rage, where we first average P(e, T) 

N = 64 typically we have Ms 10 5 . In all cases we used Qver the disorder ( aver age denoted by (T)), and thc 

N mn = 10 6 . In general, thc study of larger sizes requires configurational entropy per spin is estimated as 

a very large amount of memory and CPU time though it 

, i,, ■„ ,. ■ ! -1 • , s c (e) ee ° g ' ^' U. +p e + const, (28) 

should be possible, as we will discuss m more detail m the N 
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Fig. 1. The complexity as a function of the IS energy for the 
SK model with N = 64 at temperatures T = 1.0,0.8,0.6, for 
Gaussian and ±1 couplings. The filled symbols are quenched 
averages, the empty ones correspond to annealed averages. The 
dashed line is the analytic annealed result |29j . 



the quenched average, where 



s c (e) = l0g ^ e ' r)) + (3e + const. (29) 



In either case we have determined the normalization 
constant by superimposing the data to the Bray and Moore 
|29| analytical result, obtained within the annealed ap- 
proximation. A reasonable estimate of the size of the er- 
ror induced from sample to sample fluctuations is given 
from the dispersion of the curves at similar values of the 
temperature. 

The annealed average corresponds to what is known as 
white average in the calculation of the number of TAP so- 
lutions. For the SK model, a zero-temperature calculation 
shows [2HJ that the white average gives incorrect results 
below a critical energy e c = —0.672, where in fact the 



quenched average is needed. This can be well appreciated 
in figure^ Two results emerge from this figure: 

— the two types of average agree above the predicted 
value for e c , and clearly differ below. The strong dif- 
ference found between the two kind of average reflects 
the fact that the IS minima are uncorrelated only 
for e > e c whereas below the breaking of replica sym- 
metry should be taken into account. 

— above e c it is impossible to obtain a good collapse of 
the data, and they do not fit the theoretical prediction 

This effect, which is absent in the ROM [22] where 
all data collapse pretty well in a single curve, suggests 
that a complexity defined in terms of the usual SW 
decomposition does not give a good description of the 
free energy landscape of the SK model. 

To compute the complexity as function of IS free en- 
ergy we have estimated the IS free energy from ff 1 51) and 
evaluated the complexity using (|13fl . We note that the fi- 
nite number of searches introduces some ambiguity in the 
normalization constant. Indeed since N mn is finite the IS 
with probability 

1 



Pis CO < rf s = exp(-/3 N(f» s - f cq (T)) 



N Tl 



(30) 



i.e., with free energy /is > / IS , are never found. Fur- 
thermore, one can assume that pis(T) is correctly eval- 
uated only for IS that have been found at least few times 
(say 5 times). Consequently P(f,T) is only known for 
/ < /is anc ^ this introduces an unknown constant in s c . 
To eliminate this ambiguity we have fixed the unknown 
constant for each sample by the use of l|24|). i.e. by impos- 
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0.25 



0.15 




0.05 



Fig. 2. Complexity for the SK model as function of the free en- 
ergy for T=1.0 (circles), 0.8 (squares), 0.6 (diamonds). Empty 
symbols correspond to N = 32, filled symbols to N = 64. The 
dashed line is the analytic result. [29] In the inset we plot the 
potential <5(/, T) — f eq (T) as function of / — f eq . 

ing f e q(T) = min/#(/,T). In all cases the correction to 
s c {f, T) obtained directly from (|13[l was found very small, 
giving us further confidence in the quality of our sampling. 
Finally we note that in this case no remarkable difference 
between the annealed and the quenched average has been 
observed. 

In Figures |21 and |21 we show the complexity s c (f,T) as 
a function of the free energy for both the ROM and the 
SK model, compared with the corresponding s c (e), i.e. the 
T — > limit computed analytically in Ref. [2H] and |3U| . 
respectively. In the SK model the complexity strongly de- 
pends on T. On the other hand for the ROM it is nearly 
T-independent and, moreover, curves at different T are re- 
markably parallel each other and to the zero-temperature 
limit, in agreement with eq. (]21[l . This behavior confirms 



0.25 



0.2 - 



0.15 - 



0.1 



0.05 




Fig. 3. Same as in figureHfor the ROM. Here T = 0.8 (circles), 
0.7 (squares), 0.6 (diamonds), 0.5 (triangles), and N = 64. The 
dashed line is the analytic result |30|. 

that for the ROM the intra-valley entropy is very small, 
giving further justification to the results obtained in [22] 
using the SW decomposition. 

We note that finite-size corrections are small for s c (f, T) 
as shown by the similarity of the results for N = 32 and 
N = 64. This similarity is at variance with the strong 
finite-size effects usually found for other quantities. For 
instance one finds that finite size corrections to the P(q) 
in magnetic field in the SK model are enough strong to 
make the two peak structure in the glassy phase hardly 
visible for sizes as large as N = 1024 |3J (see also Ref. 
\&2\). However one should note that we are considering 
probability distributions of valley dependent quantities so 
that finite-size effects directly affect the normalization fac- 
tor in #(/, T) (that results in an overall constant shift in 
s c (/, T)). This effect is further reduced by fixing the con- 
stant in the complexity self-consistently through eq. I|24[) . 
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It seems therefore that the method could be usefully ap- scales are strongly separated. It moreover gives reasonable 

plied also to other system such as structural glasses, where evaluations of the relevant dynamic and static transition 

N > 60 particles are known to be enough for a reasonable temperatures. 

evaluation of s c (for a detailed study of finite size effects In summaryj we have proposed a new method to corn- 
in a glass forming model see Ref. [H]). putc tlle con figurational entropy as a function of the free 
In each figure the inset shows the (averaged) shape energy and tem perature. It provides a way to investigate 
for the potential #(/, T) at different temperatures. As is the free energy landscape of glassy systems which can be 
clearly seen the shape is different for the two models. For generally applied whenever there is a flat measure describ- 
the SK model *(/, T) has a minimum f*(T) close to the ing the prob ability to dynamically explore free energy val- 
lower (measured) bound over the free energy support in leys ^ It can be used to determ i n e the type of glassy tran- 
the whole T-range considered. The value of f*(T) con- gition (by comput ing the shape of the potential the 
verges (for N large) to / eq (T) corresponding to s c (/*) = 0. mode coupling temperature and the Kauzmann tempera- 
The ROM shows a different behavior. Here a local mini- tures 
mum appears very close to the mode-coupling transition 

The main advantage of method is that it allows to 

temperature Tmct = 0.6, and moves to lower free en- .... . . . 

evaluate the complexity without using any approximation 

ergies as the temperature is lowered |20| . The minimum 



f*(T) sticks to f e q(T) at the Kauzmann or static transi- 



for the intra- valley entropy. Moreover, one directly obtains 
an evaluation of the generalized free energy <P which con- 



tion, where s c (/*(T)) = 0. From the inset of figure El a . 

tarns the physics of the model. Un the other hand, one 

simple linear extrapolation to zero of f*(T) — f c „(T) as 

should note that the method is practically limited to rela- 

function of T gives the estimate Tr- ~ 0.26 in excellent . . . 

tive small size systems. Nevertheless, it is general enough 

agreement with the theoretical value 0.25 RU1- . . , . , , , , 

to be applied to other models of structural glasses and 

spin glasses. It would be particularly interesting to con- 

6 Conclusions sider short-range spin-glass models as well as Lennard- 

The present FEIS approach, i.e. the SW decomposition Jones § lasses ' to achieve a better understanding of the 

based on inherent structures but considered in terms of free ener gy landscape of generic glassy systems, 

their free energies, do capture the physics of the system In particular in structural glasses the total number of 

even for finite T and allows for distinguish between con- IS is still given by l|27|) with a depending on the particular 

tinuous and one step RSB scenarios. It seems particularly model and on the density, but one finds 0.7 < a < 1.2 for 

intriguing to us that it describes well also FRSB models, a Lennard- Jones model with N = 60 JS], which gives 

where it does not correspond to a situation where time- a very large IS number. Nevertheless it should be still 
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possible to carry out the present kind of analysis by fixing 3. F. H. Stillinger and T. A. Weber, Phys. Rev. A 25, (1982) 

an appropriate cutoff on pjs, which simply means that 978; Science 225, (1984) 983. 

one is estimating s c (f) for / < /° s instead of on the whole 4 - Th - M - Nieuwenhuizen, J. Phys. A 31, (1998) L201; G. Biroli 

free energy range and ^ Kurchan, Phys. Rev. E 64, (2001) 16101. 

As a last remark, it is interesting to note that at vari- 5 ' R ' G ' Palmer ' Adv - 31 ' ( 1982 ) 669 ' 

.., .. ,, ,. T „ . .... 6. J. P. Bouchaud at al., in the second reference in ffl- 
ance with the present result the lb decomposition is not 

relevant for describing coarsening models |33|. since sys- 
tems sharing the same s c (e) display a different dynamical 
behavior. 
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